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Quantification by Inhibition of Amplification 

5 The present invention relates to a method of quantifying nucleic acid 
sequences differing in at least one base using Real Time PCR. 

Background of the Invention 
Screening for point mutations, small base pair deletions or insertions in 

10 genes causing diseases has lately gained growing interest especially 
considering the development of certain diseases among an ethnic group or 
the predisposition of an individual patient. A variety of mutation scanning 
techniques (Beaudet et al, 1995, "Genetics, biochemistry and molecular basis 
of variant human phenotypes", The Metabolic and Molecular Bases of 

15 Inherited Disease, Vol. 1, McGraw-Hill, Inc., New York) with, mutation 
detection rates up to and above 90% (Guo et aU 1997, "Enhanced 
discrimination of single nucleotide polymorphisms by artificial mismatch 
hybridisation". Nature Biotechnology, Vol. 15: 331-335) have been 
developed over the last decade. 

20 However, even more of interest than the knowledge about the presence of a 
specific mutation in a nucleic acid sequence would be the knowledge about 
the amount of such mutated nucleic acid sequences. 

The standard technique to quantify nucleic acid sequences with mutations is 
a notoriously laborious combination of a polymerase chain reaction (PCR), 

25 a labeled (e.g. radioactive) restriction enzyme digestion and phospho- 
imager quantification (Saperstein & Nickerson, 1991, Biotechniques 10 (4) 
488-489). Accordingly, the sequences to quantify are amplified in a PCR, 
the amplified products are incubated with restriction enzymes and 
radioactive labeled. Only the nucleic acid sequences that carry a specific 

30 mutation, thereby comprising an additional restriction enzyme recognition 
site, can be cleaved with the restriction enzyme. Such a restriction enzyme 
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digest results in uncleaved amplification products, not comprising the 
mutation, and cleaved amplification products, because of the mutation. 
Depending on their size these products are separated by gel-electrophoresis 
separation. This gel or a membrane after southern blotting this gel is 
5 exposed to screen sensitive to radioactivity. Since all fragments are equally 
labeled by radioactivity the stimulation, which correlates with the amount 
of amplification product and which is caused by the radioactive 
disintegration, can be detected. Thus, it is possible to quantify said 
stimulation and the amount of the different amplification products, 

10 respectively, with a phospho-imager. 

Beside the time consuming setup the limits of the detection range for 
quantification are reached if less than 5% of the allelic sequences of the 
genome carry a mutation. In other words, if only 1% or less of the nucleic 
acid sequences carries a mutation, this mutation is not detectable and, thus, 

15 can not be quantified. 

Additionally, not every mutation can be detected by this technique, because 
restriction enzymes that are able to distinguish between sequences with or 
without mutation are often not available. 

20 

Recently, considering the fact that PCR technology has been proven to be a 
powerful tool for nucleic acid analysis, a highly elaborate quantification 
assay for PCR products, the so called TaqMan"^^ PCR (PE Applied 
Biosystems, Germany), was suggested (Livak et al, 1995, "Oligonucleotides 
25 with fluorescent dyes at opposite ends providing a quenched probe system 
useful for detecting PCR product and nucleic acid hybridisation", PCR 
Methods Applic. 4: 357-362). 

The TaqMan'^'^-technology is based on PCR wherein an additional specific, 
30 internal, fluorogenically labeled oligonucleotide probe is used. Typical 
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amplicons range in size from 100 bp to 1200 bp. The probe specifically 
anneals between the forward (5') and reverse (3') PCR primer binding sites 
(Figure 1). It consists of an oligonucleotide with a 5'-reporter dye (e.g. FAM, 
6<arboxyfluorescein) and a quencher dye (e.g. TAMRA, 6- 
5 carboxytetramethylrhodamine) which compensates the emission spectra of 
the reporter dye as long as both dyes are attached to the probe. While 
performing the standard PCR the 5'->3' exonuclease activity of 
conventional thermostable Taq polymerase is exploited, which degrades the 
internal fluorogenic probe thereby releasing the fluorogenic signal of the 

10 5'-reporter dye. Thus, the fluorescence signals are detectable and can be 
quantified. An ABI PRISMT^ Sequence Detection System (PE Applied 
Biosystems, Germany) is normally used to detect the fluorescent signals. 
An improvement of the TaqMan^-technoIogy was the introduction of, 
instead of an endpoint measurement, the measurement of the released 

15 fluorescent emission continuously during the PCR amplification (Heid et 
al, 1996, "Real Time Quantitative PCR", Genome Research, Vol. 6, 986-994). 
Since the exponential accumulation of the fluorescent signal directly reflects 
the exponential accumulation of the PCR amplification product this 
reaction is monitored in real time. From the output data of a so-called Real 

20 Time PCR quantification a reliable back calculation to the input target DNA 
sequence is possible. 

For the detection of minimal alterations on nucleic acid sequences found 
within alleles, genetic polymorphisms or the resulting transcripts the 

25 TaqMan'f^ technology can be adjusted. In this case, for every allele or 
polymorphism to be detected by TaqMan™ technology a specific 
oligonucleotide probe must be designed and used. These oligonucleotide 
probes differ in one or more nucleotides in accordance with the allele the 
probe shall bind to for TaqMan'^^ PCR assay. Additionally, for every allele 

30 the corresponding oligonucleotide probes comprise a different colored 
fluorescent reporter dye. Thus, by performing the TaqMan^w ^gsay 
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every allele or polymorphism is recognized by the probe that specifically 
binds to its sequence. During PCR, the fluorogenic signal of the specific 
probe is released. Since the released fluorogenic signals are different for 
every allele it can be determined how many different alleles are present. 
5 However, to determine the quantity of each different allele this method is 
unreliable or even not applicable: 

Assuming, that alleles only differ in one base the corresponding 
fluorogenic probes used for TaqMan'^^ PCR assay also differ in only one 
base. In this case the specificity of the quantification assay is due to an 

10 inhibition of hybridization. Accordingly, the fluorogenic probes are of 
minor specificity: Since a probe specific for one allele differs from another 
allele only in one base this minor difference will not prevent binding of 
said probe to the other - unspecific - allele. Thus, it happens that in a 
TaqMan™ PCR assay the same fluorogenic signal will be released from 

15 different alleles. Accordingly, this fluorogenic signal can not be used for a 
reliable quantification of the corresponding allele. 

To summarize technologies as known from the prior art are not suitable to 
quantify alleles or nucleic acid sequences with point mutations, small base 
20 pair deletions or insertions, especially, when said sequences represent less 
than 10% of the allelic sequences in the genome. 

Object of the invention 

25 

It is, thus, an object of the present invention to provide a highly sensitive 
and reliable assay capable to distinguish and especially to quantify nucleic 
acid sequences that differ in at least one base. 

30 



The Invention 
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In its broadest aspect the present invention provides a method of 
quantifying nucleic acid sequences differing in at least one base wherein an 
allele-specific primer, forming a match with one allele but a mismatch with 

5 the other(s), is used for a Real Time PGR. The primer-binding site of the 
allele-specific primer according to the present invention is selected in a way 
that the allele-specific primer binds at the location of the sequence variation 
of the allele to be quantified. Accordingly, said primer forms a match with 
the one allele but a mismatch with the other allele(s) at the location of the 

10 sequence variation. As a result of the match said allele-specific primer 
forms with the one allele of the template nucleic acid sequences, this allele 
is amplified in a Real Time PGR assay. During the extension phase of the 
Real Time PGR assay the 5'-Exonuclease activity of the Taq polymerase 
cleaves the fluorescent reporter dye (e.g. FAM, 6-carboxyfiuorescein) from a 

15 hybridization probe and thereby releases fluorescent emission of the FAM 
at 518 nm, which can be determined. 

The mismatch can be located at any position of the primer and is normally 
due to false base pairing between the nucleotides Adenine/ Tymidin and/ or 

20 Cy tosine/ Guanidine of the primer sequence with the allelic sequence. As a 
consequence of said mismatch primer-binding stability is reduced. Due to 
said reduced binding stability of the allele-specific primer this allele is not 
amplified in the Real Time PGR assay. As long as an allele is not amplified 
no 5'-Exonuclease degradation of the labeled probe occurs. Consequently 

25 the emission spectra of the reporter dye of the hybridization probe specific 
for this allele remains quenched by the second dye (e.g. TAMRA, 6- 
carboxytetramethylrhodamine) also attached to the hybridization probe. 
Therefore, no fluorescent signal of said non-amplified allele can not be 
detected. 

30 

Nevertheless, also said non-amplified allele can be quantified. 



wo 99/40219 PCT/EP99/00648 

6 

In the case that only two alleles are to be quantified, at first the fluorescent 
emission of the amplified allele is measured, as indicated above. 
Subsequently, a further Real Time PGR is performed using non allele- 
specific primers. Thus, fluorescent signals, which are released from both 
5 alleles, can be measured. Finally, the amount of the non-amplified allele is 
calculated by subtracting the determined emission value representing the 
amount of the amplified allele from the determined emission value 
representing the amount of both alleles. The main advantage of said 
amplification of both alleles is that an internal reference value is used. 

10 

Similarly, further alleles, i.e. more than two alleles can be quantified. 
For this, the amount of the first allele as well as the amount of all alleles is 
determined in a Real Time PGR, as described above. Additionally, the 
amount of all further alleles, except one non-amplified allele which will be 

15 calculated, have to be separately determined in a Real Time PGR. For every 
allele, to be measured, a different allele-specific primer is used. These 
primers form a match only with the allele they are specific for, but form 
mismatches with all the other alleles. Irrespective of the number of alleles 
present in a sample, in a Real Time PGR assay said primers allow only 

20 amplification and, respectively, determination of an emission value of the 
allele they are specific for. Finally, the amount of the non-amplified allele 
is calculated by subtracting the separately determined emission values 
representing the amount of the first and further amplified alleles from the 
determined emission value representing the amount of all alleles. 

25 

It has been found that the quantification assay is highly sensitive in 
comparison with prior art quantification techniques. It has been 
surprisingly found that an allele, which represents only 0.01% of the allelic 
sequences in a genome, can still become amplified and quantified, 
30 respectively. Thus, the method according to the present invention for the 
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first time provides a quantification assay, which enables to quantify an 
allele that represent less than 5% of the allelic sequences in a genome. 

A further advantage of the method according to the present invention is the 
5 provision of high sample throughput obtainable by automatisation. 

Moreover, compared with the state of the art methods the method 

according to the invention avoids involvement of working steps difficult to 

standardize, such as electrophoresis for the analysis of PCR products. 

Furthermore, due to the use of allele-specific primers only a single labeled 
10 probe, instead of labeled probes for every allele, is needed. This results in 

tremendous cost reduction compared with costs arising from standard 

technology assays. 

The term ^'allele" as used, therein after, is not limited to the understanding 
15 of the variants of a gene in a diploid DNA genome (Rompp Lexikon 
Biotechnologie, Georg Thieme Verlag Stuttgart - New York, 1992), but is 
used in a broader manner. The term "allele" additionally comprises other 
nucleotide acid sequences with point mutations, deletions and/or 
insertions of one or more base pairs, such as mitochondrial DNA (mtDNA), 
20 messenger RNA (mRNA), viral DNA or RNA genomes or DNA of 
microorganisms. Consequently, the method of the present invention is 
applicable for the quantification of nucleic acid sequences of other systems. 

For example, it is possible to analyze quantitative differences in parental, 
25 heteroplasmatic or mutated variants of mtDNA using the method according 
to the present invention. The distribution of parental mtDNA is of special 
interest for pedigree analysis of domestic animals. Normally 
predominantly maternal mtDNA is present in descendants, except for 
descendants generated by in vitro fertilization or cloning techniques. For the 
30 analysis of the mtDNA composition in the cloned offspring of livestock 
species a reliable quantification for mtDNA of the progenitors is essential. 
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Further, in the context of mitochondrial diseases, such as optic neuropathy 
(LHON) and others (for review see Johns (1996) The other human genome: 
Mitochondrial DNA and disease. Nature Medicine 2: 1065-1068) a reliable 
assay to determine quantitative differences of mutated mtDN A is desirable. 
5 When considering diseases due to mtDNA defects it is essential to consider 
the phenomenon known as heteroplasmy. This describes the presence of 
both mutant and wild-type mtDNA in the same tissue, cell and possibly 
within the mitochondrion. Furthermore, the proportion of mutant to wild- 
type genome is seen to vary between tissues and also throughout hfe. 

10 Therefore, much interest has been focused on the hypothesis that 
accumulating damage to mtDNA is of importance in aging and common 
neurogenerative diseases, e. g. Parkinson's disease and Alzheimer's disease 
(Egensperger et al, 1997, "Association of the mitochondrial tRNA (A4336G) 
mutation with Alzheimer's and Parkinson's disease", Neuropathol-Appl- 

15 Neurobiol 23(4): 315-321). The method according to the invention provides 
for the first time a possibility for the determination and quantification of 
organ- and age-related levels of mutated mtDNA. As a result of such a 
quantitative determination of the mutated mtDNA therapeutic approaches 
to cure e.g. an organ-related mitochondrial disease can be applied directly 

20 to the involved organ. 

As a further example, it is also possible to quantify the different RNA 
transcripts of allelic genes. Said quantitative analysis is especially of 
interest in case that the transcripts of the different alleles, because of a 

25 sequence variation on the one allele, have different functional properties. 
The term "different functional properties" comprises that an RNA transcript, 
including but not limited to mRNA, rRNA or tRNA, lost its capacity to 
induce or prevent a certain function. In this case an unbalanced or even 
predominant presence of said RNA transcript can correlate with a disease 

30 syndrome. One example for a disease caused by different functional 
properties of allelic transcripts is sickle cell-thalassemia disease which 
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depends from the percentage of the mutated Hb (Roche Lexikon Medizin, 
Eds. Hoffmann-La Roche AG und Urban & Schwarzenberg, Munchen, 1984). 
Reliable quantification of such RNA transcripts with different functional 
properties will influence decisions for further medical treatment. For 
5 quantification according to the method of the present invention, at first, 
said transcripts of allelic genes have to be reverse transcribed to cDNA. 
Subsequently an allele-specific primer designed to form a match with the 
cDNA of a certain transcript is used in a Real Time PGR quantification assay 
to amplify and, respectively, quantify the amount of said certain transcript. 
10 As already described, the amount of both transcripts, independently their 
functional differences, is determined in a Real Time PGR using non allele- 
specific primers. Finally the amount of the other non amplified transcript 
can be calculated as mentioned above. 

The method according to the present invention is also applicable for the 
determination of differences between the wildtype genome of a virus and 
mutated genomes of subtypes thereof. For example, in case of Human 
Immundeficiency Virus (HIV) it is known that resistance against a special 
drug is due to a specific mutation within the viral genome. According to the 
method of the present invention the DNA or, respectively, cDNA of the 
mutated viral genome from a patient derived sample is amplified within a 
Real Time PGR using an allele-specific primer forming a match with said 
mutated genome but a mismatch with the wildtype genome, i.e. those 
viruses not comprising the specific mutation. After an additional Real Time 
PGR using non allele-specific primers to determine the amount of all 
viruses in the patient derived sample, the amount of the non amplified and, 
thus, not mutated viral genomes can be calculated as described above. This 
quantification of the different viral subtypes in the virus population of a 
patient helps to decide on a drug therapy, which said virus population is 
still sensitive to. 



15 



20 



25 



30 
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Furthermore, the method of the present invention is applicable to quantify 
alterations in DNA of microorganisms, including extrachromosomal DNA 
of bacteria, cyanobacteria or fungi. Extrachromosomal DNA of 
microorganisms regularly comprises sequences encoding for proteins 

5 leading to a drug resistance or also toxic proteins. One example for a 
bacterium encoding a toxic protein on extrachromosomal DNA is E. coli 
subtype EHEC that can contaminate dairy products, raw meet or similar, 
and causes severe diarrhea or renal failure which leads quite often to death 
of the patient. Therefore, it is of special interest for the food industry to 

10 quantify the percentage of bacteria harboring DNA encoding for toxin 
genes. According to the method of the present invention the amount of 
bacteria carrying extrachromosomal DNA encoding for toxin genes in a 
bacterial population can be quantified. For this, the extrachromosomal 
DNA encoding said toxin gene is amplified within a Real Time PCR using 

15 an allele-specific primer forming a match with said gene but a mismatch 
with corresponding sequences not encoding said gene. After a further Real 
Time PCR using non allele-specific primers to quantify the amount of 
extrachromosomal DNA present in the sample, the percentage of 
extrachromosomal DNA encoding for said toxin genes can be calculated as 

20 described above. Since the amount of extrachromosomal DNA per 
bacterium varies a further Real Time PCR quantifying a single copy gene 
encoded on the extrachromosomal DNA is performed. In comparison with 
this further Real Time PCR quantifying a single copy gene encoded on the 
extrachromosomal DNA the percentage of bacteria encoding said toxin 

25 gene is determined. Similarly, the method according the present invention 
can be applied to quantify the percentage of bacteria harboring 
extrachromosomal DNA encoding for a protein leading to a drug 
resistance. Especially for a proper treatment of a bacterial infection with 
antibiotics it is important to know whether a patients bacterial population 

30 is resistant to a certain drug. For this, the extrachromosomal DNA encoding 
said resistance gene is amplified within a Real Time PCR using an allele- 
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specific primer forming a match with said gene but a mismatch with 
corresponding sequences not encoding said gene. In a second Real Time 
PCR non allele-specific primers are used to quantify the amount of 
extrachromosomal DNA present in the sample. The percentage of bacteria 
5 encoding for said resistance gene can be calculated as described above. As a 
result of this quantitative analysis of the amount of bacteria resistant to a 
certain antibiotic treatment can be adjusted. Another example of 
microorganism encoding for a certain toxin gene is Pastorella, especially P. 
multocida. A Pastorella infection with a patogenic strain of Pastorella e.g. in 

10 even-toed Ungulates leads to haemoragies and within 3 to 4 days finally to 
death of the animal. Therefore, there is a high interest to recognize strain 
specific differences of Pastorella. According to the method of the present 
invention a quantification of such strain specific differences is possible 
using an allele-specific primer specific to the pathogenic strain. As 

15 described before, the amount of the whole population of microorganisms is 
determined in a Real Time PCR using conventional primers. Additionally, 
the amount of a specific strain of Pastorella is quantified in a Real Time PCR 
using allele-specific primers. 

20 In a preferred embodiment of the invention the primer is designed to form 
a mismatch at the terminal 3' nucleotide (Figure 1) with the allele not to be 
amplified. This mismatch leads to an instability of the 3' end of the primer, 
resulting in a not fully hybridized 3' end. Consequently, this primer 
forming said mismatch is not elongated in a PCR and, respectively, the 

25 allele is not amplified in a Real Time PCR. 

To further ensure that no amplification of the allele that form a mismatch 
with the primer occurs, at least one additional mismatch is introduced into 
the sequence of the allele-specific primer. The term "additional mismatch" 
30 comprises single and multiple base pair exchanges, insertions or deletions 
as well as chemical modifications of bases to the primer sequence. Such 
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modifications disturb the natural occurring base pairing. Consequently, the 
introduction of such modifications e.g. introduction of a base analogue such 
as 3-nitropyrrole, a derivative of the nucleotide such as 
Xanthine/ Hypoxanthine or derivative of the nucleoside such as Inosine 
5 leads to an additional mismatch of the allele-specific primer with both 
alleles. The position of said additional mismatch can be located everywhere 
in the primer sequence. An additional mismatch leads to enhance primer 
instability. Due to this increased primer instability no amplification of said 
allele takes place in the Real Time PCR assay according to the invention. 
10 However, said additional mismatch also weakens the binding between 
primer and the allele, which forms a match with said primer and which 
should be amplified. It has, nevertheless, been found that Real Time PCR 
amplification still occurs. 

15 According to a further embodiment of the present invention, the additional 
mismatch is positioned between nucleotide -2 and -10 of the 3' end of the 
primer. Preferably, the additional mismatch is a single base pair mismatch 
at nucleotide -3 of the 3' end of the primer. Alternatively, primers with two 
or more mismatches between position -2 to -10 of the 3' end can be 

20 designed. 

The allele-specific primer can form the mismatch(es) with the major allele 
of the population. Usually, the major allele is defined as the wild type or 
naturally occurring allele, which normally is predominantly present in a 
25 sample compared to the whole population of alleles. The minor allele is 
defined as the allele less present in said sample^ which e.g. comprises an 
alteration or mutation. In another further embodiment according to the 
invention the allele-specific primer forms a match with said minor allele, 
which subsequently is amplified and determined. 



30 
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The allele-specific primers (Table 1) ASS (SeqID.: 1), AS6 (SeqlD.: 2), AS7 
(SeqID.: 3).PIRA (SeqID,: 4), ASl (SeqlD.: 11), AS2 (SeqID.: 12), AS3 (SeqID.: 
13) or AS4 (SeqID.: 14) are especially used according to the method of the 
present invention for quantification of alleles comprising mitochondrial 

5 DNA (mtDNA) sequences differing in at least one base. Additionally, for 
quantification of alleles of the mtDNA the labeled oligonucleotide probe 
IMP (SeqID.: 4) is used. Conventional primers TMl (SeqID.: 9), TM2 (SeqID.: 
10), COl (SeqID.: 6), C02 (SeqID.: 7) or C03 (SeqID.: 8) are especially used 
for the quantification of the whole amount of alleles, which resembles as a 

10 value of reference 100% of alleles present in a specific system. 

Summary of the Invention 

15 To achieve the foregoing and other objects, the present invention inter alia 
comprises the following, alone or in combination: 

A method of quantifying alleles comprising nucleic acid sequences 
differing in at least one base wherein an allele-specific primer, forming a 
20 match with one allele but a mismatch with the other(s), is used for a real 
time PGR resulting in amplification of the allele that forms a match with 
the primer followed by quantification of said amplified allele as well as of 
the non-amplified allele; 

25 the method as above, wherein said non-amplified allele forms a mismatch 
with the terminal 3' nucleotide of the allele-specific primer; 

the method as above, wherein said allele specific primer comprises 
additionally at least one mismatch; 

30 
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the method as above, wherein said additional mismatch is at position -2 to - 
10 of the 3' tail of said allele specific primer; 

the method as any above, wherein said allele specific-primer forms a match 
5 with a minor allele present in smaller amounts than the non-amplified 
allele(s); 

allele-specific primers with the sequence 





ASS 


(SeqID.:!); 


10 


AS6 


(SeqID.: 2); 




AS7 


(SeqID.: 3); 




PIRAl 


(SeqID.: 5), 




ASl 


(SeqID.: 11), 




AS2 


(SeqID.: 12), 


15 


AS3 


(SeqID.: 13), or 




ASi 


(SeqID.: 14); 



allele-specific primers as above used for a method of quantification of 
alleles comprising mitochondrial DNA sequences differing in at least one 
20 base. 

The following figures and examples will further illustrate the present 
invention. It will be well understood by a person skilled in the art that the 
25 provided figures or examples in no way may be interpreted in a way that 
limits the applicability of the technology provided by the present invention 
to this examples. 

Figure 1 is a schematic presentation of two alleles of the bovine 
30 mtDNA, which differ in one nucleotide in position 16057. The nucleotide 
positions of the mtDNA, herein after, refers to Anderson et al. (1982, 
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"Complete sequence of bovine mitochondrial DNA conserved features of 
the mammalian mitochondrial genome", J. Mol. Biol. 156: 683-717). The 
primer binding site of the forward (1) and reverse (2) primer are indicated. 
Stars mark the localization of sequence dependent mismatches between the 
5 primer and the template. Further indicated is the binding region of the 
oligonucleotide probe (5) with its attached fluorescent reporter dye (3) and 
quenching dye (4). 

Example I 

10 

Cytoplast-blastomere fusion and related embryo cloning techniques are 
used for the production of elite embryos in livestock species. Thus, the 
parental mitochondrial DNA (mtDNA) in early cloned embryos (n=26) and 
its transmission and segregation in born animals obtained by infmspecific 
15 cytoplast-blastomere fusion in Bos taurus was investigated. 

During oogenesis the mtDNA copy number is reduced to nearly one per 
organelle. Mitochondrial DNA replication will repopulate each organelle, 
sometime after the blastocyst stage and lead to relatively pure populations 

20 of each genotype that pre-existed in the original parental organelle. 

Cloning by cytoplast-blastomere fusion is based on the transfer of whole 
blastomeres, which donor the nucleus (donor blastomere = DB) from 
cleavage stage embryos to enucleated meiotic metaphase II oocytes 
(recipient cytoplast = RC). The cytoplast which includes cytoplasmic 

25 organelles like RER, Golgi, SER, tubules and other maternal factors being 
important for cell-cell communication and oocyte development is involved 
in the reprogramming of the introduced nucleus. The cytoplast's mtDNA 
type is assumed to be transmitted to the cloned offspring (Plante, 1992, 
"Restriction fragment length polymorphism in the mitochondrial DNA of 

30 cloned cattle", Theriogenology 38: 897-904). However, in the cloning 
process there are two sources, i. e. cytoplast- and blastomere-, of parental 
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mtDNA in the context of one nuclear genome, which is surrounded by two 
types of cytoplasm. 

To investigate the transmission of the mtDNA allelic differences in the 
5 bovine mtDNA control region was quantified. 

At first, and to have a comparison both types of parental mtDNA in cloned 
embryos of early developmental stages was determined by Primer 
Introduced Restriction Analysis (PIRA). The A to G transition at ntl6022 
allowed the differentiation of parental mtDNA used (nucleotide positions 

10 refer to Anderson et al, 1982). The rare allele 16022A present in the DB- 
mtDNA created a potential restriction site for Nspl which was used for 
primer introduced restriction analysis. The primer PIRAl (SeqID No.: 5) 
possessing a mismatch at position 3 from the 3' end, generated a Nspl site 
in the PGR product amplified with PIRAl (SeqID No.: 5) and C02 (SeqID 

15 No.: 7). PGR amplification was carried out in a buffer containing 60 mM 
Tris-HCl (pH 8.5), 15 mM (NH4)2S04, 0.2 mM dNTPs, 1 ^M of each primer, 2 
units Taq polymerase (Life Technologies, Austria) and 1 |il template DNA. 
PGR amplification was performed in a 50 ^il reaction volume on a 
RoboCycler® (Stratagene, USA). The amplification of embryonic material 

20 was carried out under hot start conditions applying Hot-wax beads® 
(Invitrogen, The Netherlands) which contained MgCh giving rise to final 
concentrations of T5 mM. Each experiment included a NTC (no template 
control). The amplification was carried out routinely according to the 
following scheme: first cycle: 5 min 95 °C/ 1 min at annealing temperature 

25 (Ta) ^VC/ 3 min at 7TC) 38 cycles: 20 sec 95°C, 40 sec at Ta/ and 60 sec 7TC 
and last cycle: 20 sec 95^C 40 sec at Ta/ 180 sec 72°C, As a control for the 
intactness of template DNA the whole control region (GR) of bovine 
mtDNA was amplified with primers COl (SeqID No.: 6) and G03 (SeqID 
No.: 8), performing a hot-start PGR with Ta ^ 48 °C 1.5 mM MgGh. 
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In the last two cycles 1 \iCi of [(x32]P-dCTP was added to label the PGR 
product. 15 1x1 of this reaction was cut with 2 U Nspl (Life Technologies, 
Austria) for 3 h and run on a 10% non-denaturing polyacrylamide gel. The 
DB-mtDNA was taken as a control for a complete restriction digest. Gels 
were analysed on a Biorad MG525 phospho-imager (Bio-Rad Laboratories, 
Austria). 

It was shown that until the blastocyst stage the amount of mtDNA in the 
cloned embryos was constant The percentage of the DB-mtDNA compared 
to the RC-mtDNA are 13 %, 13 % and 9 % for the 1-cell stage, the morula 
and the blastocyst, respectively. 

Subsequently the transmission of parental mtDNA in three pairs of cloned 
calves after birth was assayed by allele-specific quantification according to 
the present invention. Therefore, clones were generated from donor 
15 embryos with differing numbers of blastomeres: 24-cells, 52-cells, and 92- 
cells. Total DNA was isolated from blood as described by Kawasaki (1990). 
AUele-specific primers ASS (SeqID No.: 1), AS6 (SeqID.: 2) and AS7 (SeqID 
No.: 3) were designed instead of conventional oligonucleotides to exclude a 
possible amplification of a false 'mitochondrial allele'. Therefore, an 
20 additional mismatch at position 3 from the 3' end of each primer was 
introduced. 

The probe TMP (SeqID No.: 4) for the allele-specific TaqMan^w 
quantification was designed to anneal between the forward and reverse 
primer (Figure 1). It consists of an oligonucleotide with a 5' -reporter dye 
25 (FAM, 6-carboxyfluorescein) and a quencher dye (TAMRA, 6- 
carboxytetramethylrhodamine) attached via a linker arm that is typically 
located at the 3' end. The probe and the allele-specific primers were of 
HPLC-grade. A Hot Start PGR was not necessary under these conditions. 
TaqMan'f'^ PGR was performed at 4 mM MgCh (5 min at 95°G, 40 cycles for 



5 



10 
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15 sec at 95^C and 1 min at 70^q. To test the reliability of TaqMan™ PGR in 

each experiment the amplification of each sample was repeated three times 
using the same master-mix. 

A computer algorithm compares the amount of reporter dye emission with 
5 the quenching dye emission (Q) in each step and each cycle during the PGR 
amplification, generating a ARn (also called ARQ, ARQ = R/Q) value. This 
value reflects the amount of hybridization probe that has been degraded. 
The algorithm fits an exponential function to the mean ARn values of the 

last three data points of every PGR extension cycle, generating an 
10 amplification plot. A relative fluorescent emission threshold is set based on 
the baseline of the ARn during the first 10-15 cycles. The algorithm 
calculates the cycle at which each PGR amplification reaches a significant 
(i.e. usually 10 times the standard deviation of the baseline) threshold (Gt). 
Therefore, the threshold cycle (Gt), i.e. the first cycle in which a PGR 
15 product was detected, is therefore given as a mean value. It was 
demonstrated for conventional Real Time TaqMan PGR that the calculated 
Gt value is proportional to the number of target copies present in the 
sample (Heid, etal, 1997). Thus, the Gt value is a quantitative measurement 
of the copies of the target found in any sample and the real-time PGR 
20 method allows for a very large assay dynamic range (approaching 10^-fold 
starting target). 

In addition to the NTG as mentioned above a second type of control the 
NAG (no amplification control) was included. Therefore, total cellular DNA 
of cattle, which lacked the mitochondrial allele under study, was used. If the 
25 desired type of mtDNA was not available a plasmid (pTAg vector, R&D 
Systems, Germany) bearing the appropriate allele was used. 

The ratio of parental mtDNA types was analyzed by two sets (each set in 
triplicates) of TaqMan*^^ experiments, which were necessary to quantify the 
percentage of DB-mtDNA compared to the total population of mtDNA 
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molecules: (i) the determination of the total amount of mtDNA (NAC) and 
(ii) the allele-specific TaqMan'^'^ PGR of interest. 

The calculation for allele-specific quantification can be further illustrated by 
5 the' following equation: 

n 

n=l 

10 

Zf total number of molecules; n, number of different alleles, and 
Xn, number of molecules represented by each allele. 

Instead of total numbers of molecules the allelic distribution can be given 
in percent to the entirety of molecules (z - 100 %). The determination of 
15 parameters involved in this equation is described in the following for a 
two-allelic system. The value z can be determined in a conventional Real 
Time PGR. The allele-specific quantification determines one of the two 
alleles, which is usually the less frequent one. 

To summarize it could be shown for conventional 'in vitro fertilizeci' (IVF) 
20 embryos and for cloned embryos that the content of mtDNA is constant 
during early embryogenesis until the blastocyst stage suggesting an 
absence of mtDNA replication. Therefore, the ratio of parental mtDNA was 
in accordance with the estimated quantitative participation of mtDNA from 
the fusion partners. Evidence has been obtained for a low-level 
25 transmission of blastomere mtDNA (DB-mtDNA) into the cloned offspring, 
thereby generating a heteroplasmic population of mtDNA. This has 
implications for cloning livestock animals. The amount of DB-mtDNA was 
13 % and 18 % in two animals of a clone which derived from a 24-cell 
morula (donor of blastomeres) and 0.6 % and 0.4 % in two calves of clonal 
30 origin derived from a 92-cell morula (Table 2). These values are in 
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accordance with the tendency expected for neutral mtDNA segregation that 
the less cell divisions that have occurred in the donor embryo, the higher 
the amount of DB-mtDNA. We also found a strong decrease of DB-mtDNA, 
which was about three magnitudes in the third clone derived from a 52-cell 
5 morula stage. 

Example 2 

The applicability of the present invention for the quantification of 
alterations in microorganisms/ which include chromosomal and especially 

10 extrachromosomal DNA, is demonstrated in the following for the case of 
extrachromosomal DNA. As model we have used E. co/i transformed with 
different plasmids, pl69A and pl69G. These plasmids represents a DNA 
construct consisting of a cloning vector (pTAg vector, R&D Systems, 
Germany) and the whole bovine mtDNA control region (nucleotide 15747 

15 to 383; Steinborn et al, 1998, " Genetic variation in functionally important 
domains of the bovine mtDNA control region", Biochim. Biophys, Acta.). 
The plasmids varied in the way that, pl69A contained the mtDNA allele 
169 A, which carries the nucleotide A at position 169, and pl69G contained 
the other allele 169G, carrying the nucleotide G at position 169. Plasmid 

20 DNA was isolated from the bacteria culture according to a standard 
minipreparation method. Commonly, the plasmid DNA is isolated at a 
concentration of 250 ng/|il and was subsequently diluted 10^ fold. The 
applicability of the allele-specific Real Time PCR for the quantification of 
different alleles of extrachromosomal DNA was demonstrated in 

25 experiments mimicking the in vivo situation in which a population of 
molecules contains different ratios of both plasmids. For this the 
amplifications were performed on a 1:2 serial dilution of the plasmid 
(pl69A) harboring the allele 169A in the background of plasmid pl69G, 
harboring the allele 169G. 
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For a quantification of the plasmids coding for the different alleles the 
fluorescent probe TMP (SeqID.: 4) and the allele-specific primer AS7 
(SeqID.: 3) and as reverse primer TMl (SeqID.: 9) were used. As a surprising 
effect it was shown that the method according the present invention 
5 remains highly specific even using an allele-specific primer with a 
nucleotide A as a terminal 3' end. In a second Real Time PCR the two 
conventional oligonucleotides TMl (SeqID.: 9) and TM2 (SeqID.: 10) as well 
, as the fluorescent probe TMP (SeqID.: 4) were used to quantify the whole 
amount of plasmid. The expected differences in the threshold values * an 1:8 

10 dilution (2^ dilution steps) is expected to show three times later its 
threshold - were obtained for each of these serial dilutions. Thus, it was 
demonstrated in experiments mimicking the in vixv situation of mtDNA 
heteroplasmy, i.e. the occurrence of at least two mtDNA alleles in the same 
sample/ cell, that the allele-specific Real Time PCR is a reliable 

15 quantification assay. 

Furthermore, several parameters influencing the efficiency of Real Time 
PCR amplification were investigated. Such parameters are e.g. magnesium 
and salt concentrations, reaction conditions (i.e. time and temperature), size 

20 and composition of the expected amplification product, primer sequences, 
and sample purity. Due to the mismatch(es) of the allele-specific primer 
near its 3' end this type of primers are especially sensitive to sample purity. 
In preliminary experiments it was suggested that the quality of the DNA 
preparation might be a limiting factor for some allele-specific primers due 

25 to the specific primer sequence. Therefore, the issue of sample purity in 
relation to a possible inhibitory effect was analyzed for 10^ and 10^ 
dilutions of the original plasmid isolation. As a result of this experiment it 
was shown that no inhibitory effect using the above mentioned allele- 
specific primers for the allele-specific amplification was observed. 
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Table 1 

AUele-specific oligonucleotides 



Oiigonucleotide 


5'-3' sequence 


S'-3' position 


SeqID.: 


ASl 


gta ctt get tat atg cat ggT(c) gt 


16044-16022 


11 


AS2 


taa tta tat gta tta tgt acG(a) gg 


16079-16057 


12 


AS3 


cca gca taa tga taa C(g)ca 


152-169 


13 


AS4 


gag cac cag cat aat gat aaA(g) eg 


147-169 


14 


ASS 


ace att gac tgt aat gte T(g)at 


189-169 


1 


AS6 


gcc cca tgc ata taa gca G(a)gt 


16022-16042 


2 


AS7 


gca agt aca tga cct eta C(t)ac 


16037-16057 


3 


PIRAl 


atg tat ata gta cat taa att aC(t)a t 


15997-16021 


5 


COl 


cac cat caa ccc cca aag ct 


15747-15766 


6 


C02 


cct gaa gaa aga ace aga tg 


16284-16265 


7 


COS 


ttg ggt taa get aca tea ac 


383-364 


8 


TMl 


ctt aat tae cat gcc gcg tga 


16159-16179 


9 


TM2 


cca get aca ata gat get ccg 


131-111 


10 


TMP 


ttg acg gee ata get gag tec 


99-79 


4 



5 

The position refers to the numbering. in the reference sequence (Anderson et 
al, 1982). Nucleotides generating a mismatch in the template DNA are 
printed as capital letters. The corresponding nucleotide of the reference 
sequence is given in brackets. 
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Table 2 

Different ratio of mtDNAs in cloned cattle quantified by allele-specific Real 
Time PGR. 



clone 


sample 


amount of 
mtDNA* 


allele-specific 
Real Time PCR 


%ofDB- 
mtDNA 


C24 


DB-24-cell 


20.8 ±0.2 


24.0 ±0.0 


100 % 




C24-1 


23.1 ±0.2 


29.2 ±0.1 . 


13% 




C24-2 


20.3 ±0.1 


25.9 ±0.2 


18% 




NAClt 


20.3 ±0.1 


30.2 ±0.1 


0.9% 


C52 


DB-52.cell 


20.9 ± 0.3 


20.6 ±0.3 


100 % 




C52-1 


21.0 ±0.1 


34.8 ±1.1 


0.006 % 




C52-2 


21.0 ±0.1 


38.5 + 0.6 


0.0004 % 




NAC2 


21.7 ±0.3 


40.0 ±0.0 


0% 




NAC3 


21.0 ±0.1 


40.0 ±0.0 


0% 


C92 


DB-92-cell 


22.7 ± 0.1 


23.4 ±0.1 


100 % 




C92-1 


16.8 ±0.1 


24.8 ±0.2 


0.6 % 




C92-2 


18.2 ±0.1 


26.7 ±0.2 


0.4% 




NAC4 


23.0 ± 0.1 


40.0 ±0.0 


0% 



5 The total amount of mtDNA was determined by conventional TaqMan'^^PCR; 

sample NACl represents plasmid DNA; NAC2-NAC4 represents native 
bovine mtDNA. The numbers in column "amount of mtDNA" and "allele- 
specific Real-Time PCR" represent the threshold cycle, i.e. the threshold value 
for monitoring the reporter fluorescence. 
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Claims 



(1) A method of quantifying alleles comprising nucleic acid 
sequences differing in at least one base wherein an allele-specific 

5 primer, forming a match with one allele but a mismatch with the 

other(s), is used for a real time PCR resulting in amplification of the 
allele that forms a match with the primer followed by 
quantification of said amplified allele as well as of the non- 
amplified allele. 

10 

(2) The method according to claim 1, wherein said non-amplified 
allele forms a mismatch with the terminal 3' nucleotide of the 
allele-specific primer. 

15 (3) The method according to claim 1 or 2, wherein said allele 
specific primer comprises additionally at least one mismatch. 

(4) The method according to claim to 3, wherein said additional 
mismatch is at position -2 to -10 of the 3' tail of said allele specific 

20 primer. 

(5) The method according to anyone of the preceding claims, 
wherein said allele specific-primer forms a match with a minor 
allele present in smaller amounts than the non-amplified allele(s). 

25 

(6) Allele-specific primers with the sequence 
ASS (SeqID.: 1) 

AS6 (SeqID.: 2) 

AS7 (SeqID.: 3) 
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PIRAl (SeqID.: 5), 

ASl (SeqID.: 11), 

AS2 (SeqID.: 12), 

AS3 (SeqID.: 13), or 

AS4 (SeqID.: 14); 

(7) Allele-specific primers according to claim 6 used for a method 
of quantification of alleles comprising mitochondrial DNA 
sequences differing in at least one base. 
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SEQUENCE LISTING 



(!) GENERAL INFORMATION: 

(i) APPLICANT: 

(A) NAME: Bavarian Nordic AS 

(B) STREET: Naverland 2 

(C) CITY: Glostrup 

(E) COUNTRY: Denmark 

(F) POSTAL CODE (ZIP): 2600 

(ii) TITLE OF INVENTION: Allele specific quantification 
(iii) NUMBER OF SEQUENCES: 14 

(iv) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS/MS-DOS 

(D) SOFTWARE: PatentIn Release #1.0, Version #1,30 (EPO) 



(2) INFORMATION FOR SEQ ID NO: 1: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 21 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 



(xi) SEQUENCE DESCRIPTION; SEQ ID NO: 1: 
ACCATTGACT GTAATGTCTA T 21 



(2) INFORMATION FOR SEQ ID NO: 2: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 21 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
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(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 
GCCCCATGCA TATAAGCAGG T * 21 



(2) INFORMATION FOR SEQ ID NO: 3: 



(i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 21 base pairs 
{B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc ' "primer" 

(iii) HYPOTHETICAL: YES 



Ui) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 
GCAAGTACAT GACCTCTACA C 21 



(2) INFORMATION FOR SEQ ID NO; 4: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH; 21 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4: 
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TTGACGGCCA TAGCTGAGTC C 



(2) INFORMATION FOR SEQ ID NO: 5: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 25 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 
ATGTATATAG TACATTAAAT TACAT 



(2) INFORMATION FOR SEQ ID NO: 6: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 19 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 
CACATCAACC CCCAAAGCT 



(2) INFORMATION FOR SEQ ID NO: 7: 

(i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 20 base pairs 
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(B) TYPE: nucleic acid 

(C) STRANDEDNESS ; single 

(D) TOPOLOGY; linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 7: 
CCTGAAGAAA GAACCAGATG 20 



(2) INFORMATION FOR SEQ ID NO: 8: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 8: 
TTGGGTTAAG CTACATCAAC 20 



(2) INFORMATION FOR SEQ ID NO: 9: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 21 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 



(iii) HYPOTHETICAL: NO 
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(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 9; 
CTTAATTACC ATGCCGCGTG A 21 

(2) INFORMATION FOR SEQ ID NO: 10: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 21 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: Other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 10: 
CCAGCTACAA TAGATGCTCC G 21 



(2) INFORMATION FOR SEQ ID NO: 11: 

(i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 23 base pairs 
{B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 11: 
GTACTTGCTT ATATGCATGG TGT 23 

(2) INFORMATION FOR SEQ ID NO: 12: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 23 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
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(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = -primer" 

(iii) HYPOTHETICAL; YES 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 12: 
TAATTATATG TATTATGTAC GGG 



(2) INFORMATION FOR SEQ ID-NO: 13: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 18 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 13: 

CCAGCATAAT GATAACCA 



(2) INFORMATION FOR SEQ ID NO: 14: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 23 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: other nucleic acid 

(A) DESCRIPTION: /desc = "primer" 

(iii) HYPOTHETICAL: YES 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 14: 



GAGCACCAGC ATAATGATAA ACG 
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